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We study spin parity effects and the quantum propagation of solitons (Bloch walls) in quasi-one 
dimensional ferromagnets. Within a coherent state path integral approach we derive a quantum field 
theory for nonuniform spin configurations. The effective action for the soliton position is shown to 
contain a gauge potential due to the Berry phase and a damping term caused by the interaction 
between soliton and spin waves. For temperatures below the anisotropy gap this dissipation reduces 
to a pure soliton mass renormalization. The quantum dynamics of the soliton in a periodic lattice 
or pinning potential reveals remarkable consequences of the Berry phase. For half-integer spin, 
destructive interference between opposite chiralities suppresses nearest neighbor hopping. Thus the 
Brillouin zone is halved, and for small mixing of the chiralities the dispersion reveals a surprising 
dynamical correlation: Two subsequent band minima belong to different chirality states of the 
soliton. For integer spin the Berry phase is inoperative and a simple tight-binding dispersion is 
obtained. Finally it is shown that external fields can be used to interpolate continuously between 
the Bloch wall dispersions for half-integer and integer spin. 
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I. INTRODUCTION 

Quantum effects in low-dimensional magnetism are a 
fascinating subject which has attracted much interest 
over the years. A notable example are antiferromagnetic 
chains where the quantum spin (or BerryEJ) phase leads to 
remarkable parity effects. It is for integer spin S only that 
the groundstate exhibits an excitation (or Haldand3) gap 
whereas for half-odd integral S such gaps are suppressed 
by interfering Berry phases Related to this phenomenon 
is the suppression of mesoscopic stiffness fluctuations for 
S half-integral antiferromagnets, whereas such fluctua- 
tions grow with chain size for integer S (similar to uni- 
versal conductance fluctuations in mesoscopic metals) fl 

Over the last years, the rapid advances in nanos- 
tructure technologyQ have opened the door to a new 
class of magnetic systems: Small single domain particle^ 
that display striking mesoscopic quantum phenomenacl 
(MQP) such as quantum coherence, quantum tunnel- 
ing, or spin parity effects. These particles exhibit one 
or several directions of minimal anisotropy energy be- 
tween which the spins can tunnel coherently. Motivated 
by theoretical predictions for uniform ferromagnetsErt£l 
and antiferromagnetalJ several experiments at sub- 
Kelvin temperatures have either aJacpm temperature in- 
dependent-relaxation phenomenalij^tj or a well-defined 
resonancellj (in the ac susceptibilibi] which scales expo- 
nentiaUji. with the number of spinsEJ in accordance with 
theory.llil Although these observations have been criti- 
cized on the basis of dissipation models, such as the in- 
fluence of nucleax-spinsliZI, the experiments on antiferro- 
magnetic ferritinllj provide a strong indication that the 
spins indeed tunnel at low temperatures. |— . |— ■ 

In subsequent work, it has been shownll3n22l that also 
tunneling depends on the spin parity via Berry phases, 
and that the magnetization switching is allowed only if 
the total spin of the particle is integral but not otherwise. 



found in uniform antiferromag- 



Similar results p 
netic particles.l 

While such spin parity effects are sometimes related 
to Kramers' idgseperacy, in particular in single domain 
ferromagnets0'EHl3, they typically go beyond this theo- 
rem in rather unexpected ways.HEil It is notably for non- 
uniform magnets that such effects can be quite intriguing 
as we know since Haldane's work on antiferromagnets.EI 
On the other hand, there has not been much related 
study on non-uniform ferromagnets, primarily because 
their groundstate is trivial and did not seem to offer much 
room for surprises. However, this is by no means so, and 
it is one of our goals to show that ferromagnets with more 
than one magnetic domain do exhibit novel spin parity 
effects and that these effects can lead to experimental 
consequences. 

We address the issue of spin parity in the context of 
MQP, although the Berry phase effects discussed here are 
of general relevance in low-dimensional magnetism. We 
start by considering the coherent quantum propagation of 
Bloch walls in the presence of periodic pinning potentials. 
Such potentials are naturally provided by the underlying 
crystal lattice or some superlattice structure that can be 
created by peripdir. deposition of materials with different 
anisotropies.oEllEj Parenthetically we note that periodic 
pinning provides a much smaller barrier height and tun- 
neling distance than one isolated pinning center would 
do. Thus the tunnclp|S-probability will be drastically 
enhanced in this caseEJO compai:fid^tfl the more tfi^<fe 
tional scenario where experimentcJcU'tJ and theorycS l3 
focus on wall tunneling out of single pinning centers. 

In a collective coordinate description the Bloch wall 
is then seen to behave like a single degree of freedom 
moving in a periodic structured This in turn results in 
characteristic Bloch bands in reciprocal space, where the 
bandwidth is determined by the tunneling rate through 
the potential. It is now at this stage where the Berry 
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phase enters the wall dynamics via an effective gauge po- 
tential that depends on the chirality, i.e. the internal 
rotation sense of the Bloch wall. For half-integer spins 
this gauge potential induces a halving of the associated 
Brillouin zone. At the same time a remarkable dynamical 
correlation occurs: Two subsequent band minima belong 
to opposite chiralities. Thus, the chirality of the wall al- 
ternates when the system is adiabatically driven through 
the Brillouin zones by an external magnetic field. As we 
shall argue, this phenomenon can be experimentally ob- 
served if there is a finite tunneling probability between 
the chiralities, a condition which is not difficult to meet in 
real systems. Due to the topological nature of the Berry 
phase these results are independent of details such as 
shape of the soliton and the pinning potential. Thus we 
expect that band-halving and chirality correlation also 
occurs in the limit of a spin-i chain where the soliton 
width approches one lattice constant. 

Besides these spin parity effects, the bandstructure 
leads to interesting coherence:, effects in form of Bloch 
oscillations of the wall center .Ei'tiJ As a result the sample 
magnetization oscillates in response to a static magnetic 
field, a behavior which is very similar to the ac- Josephson 
effect in superconductors. 

In principle these results hold for an arbitrary number 
Na of coupled ferromagnetic chains. However, observa- 
tion of MQP becomes increasingly difficult with increas- 
ing Na since observability requires tunneling exponents 
(which grow with Na) to be of the order of Planck's con- 
stant. This necessarily limits the size of sample cross 
sections (but not their lengths) and restricts considera- 
tions to low-dimensional ferromagnets, most typically of 
quasi-one dimensional size. An important consequence 
of this reduced dimensionality is the fact that dissipa- 
tion due to spin waves has a negligible effect on the wall 
dynamics since there is an associated finite size gap (be- 
sides the anisotropy gap) in the spin wave spectrum. It 
is due to these gaps that at temperatures typically below 
100 mK the spin waves freeze out exponentially fast, and 
are thus simply irrelevant for dissipation (they only lead 
to a minor soliton mass renormalization as we shall show 
explicitly.) 

To simplify our discussion we conside r i p .the following 
the limit of large hard-axis anisotropyHo as it occurs 
for instance in an Yttrium-Iron-Garnet (YIG) sample of 
the shape shown in Fig |^. We can then eliminate the 
out of easy-plane degree of freedom and the spin model 
reduces to that of a sine-Gordon model plus a gauge 
term coming from the Berry phases. In a quantum field 
approach we introduce collective coordinates, eliminate 
the spin waves and arrive at an effective action for the 
wall position. The spin waves give rise to a non-local 
term in the action which can be cast into the well-known 
Caldeira-Leggett form at low temperatures. In this way 
we make contact with phenomenological formulations of 
dissipiatiorOL^ as extensively discussed in the context of 
MQPlZI and show that the spectral function has a gap 
due to anisotropies. While there has been a number 



of work in jjiaxipus contexts related to intrinsic soliton 
dampingtScj'L^, we believe that the novel approach pre- 
sented here is most adequate to the combined description 
of wall dynamics and Berry phases and, moreover, pro- 
vides the first complete discussion of spin wave dissipa- 
tion, particularly in the context of MQP. Finally we note 
that a brief account oLflart of the results presented here 
has been given before.oEij. 

The outline of the paper is as follows. In Sec. || we 
discuss the derivation of the spin action plus topological 
phase from the Heisenberg model. Details of this deriva- 
tion via coherent spin states together with a unified treat- 
ment of the Berry phase in different gauges are given in 
Appendix In Sec. Ill we discuss static Bloch wall 
solutions and derive the sine-Gordon action plus gauge 
term in Sec. To gain confidence in our approach we 
first consider the uniform limit and show that thia-gauge 
term reproduces the known spin parity behavior .E3 As a 
byproduct we also obtain the tunneling prefactor. In Sec. 
^ we discuss the coupling between Bloch wall and spin 
waves and show that spin wave dissipation is negligible 
at low temperatures, technical details are presented in 
Appendix M. In Sec. VI we discuss the infiuence of the 



Berry phase on the Bloch wall dyna mics in a periodic po- 
tential, first in the nearly free (Sec. VI B| ) and then in the 
tight binding limit (Sec. VIC). In both cases it is shown 
that for half integral spin the Brillouin zone is halved 
and the chirality alternates. Experimental implications 
are given in Sec. VI], where we also give results for the 
level splitting due to the tunneling between the two wall 
chiralities. Finally in Sec. VIII we discuss how the in- 
terference effects are altered by external fields. A note 
regarding the terminology: The terms soliton and Bloch 
(or domain) wall are used interchangeably to denote the 
transition region between domains in ferromagnets. 



II. MODEL 

In this section we derive a continuum field theory 
to describe the quantum dynamics of nonuniform spin 
configurations in ferromagnets. Our starting point is 
a microscopic Heisenberg spin-Hamiltonian with local 
anisotropies. The transition amplitude between two 
nonuniform spin configurations is then expressed as a 
coherent state path integral. The corresponding action 
differs from the classical micromagnetic expression by a 
total derivative. While this term does not affect the clas- 
sical equations of motion, it gives rise to quantum me- 
chanical interference effects and thus leads in a natural 
way to the quantization of micromagnetics. Several ex- 
amples of such interference effects will be discussed below 
in Sections |rv|Jvl-[VIH. 



Ferromagnetic insulators can often be described by a 
Heisenberg Hamiltonian with anisotropies 
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where denotes the spin operator at the lattice site i. 
For simpUcity we assume that the spins form a simple 
cubic lattice of lattice constant a. Throughout this work 
we shall use units such that Ti — 1. The first term on the 



rhs of (2.1) is the exchange interaction between a spin at 
the lattice site i and its nearest neighbors at the lattice 
sites i + p. The next term is an easy-axis anisotropy along 
the y-axis with anisotropy constant Ky > 0. The third 
term is a hard-axis anisotropy of strength > which 
renders the xy-plane an easy plane. The spins will thus 
preferably point parallel or antiparallel to the y-axis. 

The anisotropies that are used in (2.1) are effective 
anisotropiesCJ and may arise from two different micro- 
scopic mechanisms. One contribution is the magne- 
tocrystalline anisotropy which is due to the interaction 
of the magnetic moments with their neighboring atoms 
via spin-orbit interaction. Consequently this contribu- 
tion reflects the symmetry of the crystal lattice. The 
second contribution is the dipolar interaction between 
the magnetic moments. Due to its long range nature 
this contribution depends on the sample shape and is in 
general a nonlocal functional of the magnetization con- 
figuration. It is this magnetostatic interaction that gives 
rise to the existence of domains in macroscopic samples. 
However, for quasi-one dimensional configurations this 
interaction considerably si mpl ifies and can be modeled 
by local anisotropies as in (2.1). 

Our focus in this work will be on elongated samples 
as shown in Fig. ^ with transverse dimensions smaller 
than the length scale QlJ|/lftrz]^/^. Spin waves running 
transverse to the sampleHcJ then exhibit a finite aiae gap 
such that they are frozen out at low temperatures .c3 This 
condition is met in most experimental situations stud- 
ied so far and thus we shall use a quasi-one dimensional 
model in the following. Truly three dimensional samples 
where all degrees of freedom are allowed to be excited are 
of rather limited interest for MQP since their tunneling 
rates and associated crossover temperatures (separating 
classical from quantum regime) are in general too small 
to be observedH 

We now turn to the path-integral form ulation of the 
system described by|-the Hamiltonian ( |2.l| ). We introduce 
coherent spin stateso at each lattice site, defined by fli ■ 
Sil^ti) — S\fli) where = (sin 6* cos 0, sin 6* sin 0, cos 61) 
is a unit vector. The whole system is then described by 
a product of coherent states at each of the Nl lattice 
sites, i.e. |{f2}) — ^fj^i Since we are interested 

in configurations that are varying slowly compared to 
the lattice constant, the spin state can be described by 
a smoothly varying unit vector field ft{x,T) depending 
on the coordinate x along the sample and the imaginary 
time T. As outlined in Appendix A, the transition am- 
plitude between the two states |{r2a}), KJ^b}) can then 
be expressed as an (imaginary) time path integral, 

({Jlb}|e-^«|{f24) = y2?0P(cos%-5-[^^''l, (2.2) 



where the integration is over all configurations that sat- 
isfy the boundary conditions ^{x, 0) = fla{x), ^{x, (3) — 
Of, (a;) (spatial boundary conditions will be specified 
later). The Euclidean action is given by 



Se — S^ 
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dr H, 



(2.3) 



with the Wess-Zumino or Berry phase term 



Swz — 



.SNa 



dr 



L/2 



L/2 



dx4>{l~cosB), (2.4) 



with (j) = dr4> and where a gauge has been chosen with 
the coherent states underlying the path integral (2.2) 
expressed in the "no rth-pole" parametrization (cf. Ap- 
pendix|^. Eq. (2^) has for closed trajectories the form 
of the sum over the Berry phases of all NaL/cl spins. The 
energy is given by 



H = N, 



[ ' dx\ 

J-L/2 [ 

Ky [sin^ 6 sin^ <; 



j[{d^ef + sin2 e{d^^)^ 



Kz cos 



(2.5) 



Na denotes the number of spins in the cross sectional 
area A of the sam ple, and L is the sample length. The 
parameters in (2.5) are related to those in ( pT| ) via 



J = JS^a, Ky^z = ky^zS'^/a. 



(2.6) 



The energy (2.5) is identicaLto-the traditional micromag- 
netic energy expression.oc3cj J and Kyz can now be 
related to the micromagnetic anisotropy and exchange 
constants J = Aa? , Ky = K^a^, Kz — KhO^- For an 
elongated slab as shown in Fig. |^, we have K^ = Ke^crysu 
Kfi = ifh,cryst + 2-kMq, where Mq = g^sS/a^ is the 
saturation magnetization, and i^e.cryst, -?^h,cryst describe 
crystalline anisotropies. Note, however, that for other 
sample geometries the demagnetizing term enters in a 
different form. For instance, for a cylindrical wire we 



would have Kp 



K 



e,cryst 



ttMq, while the hard-axis 



anisotropy would be of purely crystalline origin. [For 
other examples see Fig. 8 of Ref. The demagne- 

tizing energy is not always important, in particular for 
samples with misoriented anisotropy axes [see p. 15 of 
Ref. ^ the crystalline anisotropies can be much larger 
than 27rM|^. 

In saddle point approximation, 5Se = 0, and rotating 
to real time t = —it, we recover the classical Landau- 
Lifshitz equations of motion 



sinO dfcj) 
dt9 



a SH 

~ S ~Sd' 
a 1 SH 



S sine 



(2.7) 
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These classical equa tions are not affected by the total 
derivative in p. J X and thus follow from the classical 
Lagrangian densityEa £ = —^^^dt4'Cos9 + H. 

Note, however, that the (p term is of crucial impor- 
tance for the quantum dynamics: While the path inte- 
gral (2.2) contains higher winding number contributions 



(±cx)) = 0, this can be 



where a path retraces itself, the Wess-Zumino term en- 
forces quantization by destructive interference of paths 
which do not satisfy the condition S'X^i ^« = 27rn, where 
Ai is the area enclosed by the trajectory of the i-th spin 
on the unit sphere. If the (j) term were dropped in (2.4) 

— a "gauge" that has sometimes been used in the lit- 
erature — the area Ai would be measured with respect 
to the equator and one would have to impose the addi- 
tional constraintE^I that the paths not intersect the "date- 
line" . This constraint is very difficult to handle within 
a path integral formalism. On the other hand, ignoring 
this constraint would lead to a wrong semiclassical quan- 
tization of half-integral spins. Moreover, one would not 
obtain the suppression of tunnelixig for half-integer spins 
in small ferromagnetic particlestj, in clear contradiction 
to Kramers' theorem which requires that the ground state 
not be split. In Appendix A we show that all these 
difficulties can be avoided if one starts from one single 
premise — the single-valuedness of the coherent states 

— which leads to a restricted set of "admissible" gauges. 
Finally, we remark that if we work in the south-pole 

parametr izat ion of the coherent state (cf. |A^ ), the <j)- 
term in (2.4) changes sign but, of course, all physical 
effects that will be derived below are independent of the 
gauge (provided the gauge is admissible). 



III. BLOCK WALL CONFIGURATIONS 

There are two energetically dege ner ate spin configura- 
tions which minimize the energy (2.1): Uniform config- 
urations with all spins pointing either along the positive 
or along the negative j;-direction. We are now interested 
in structures that interpolate between these two config- 
urations. Due to the easy-axis anisotropy in (2.1), this 
transition region will have a finite width and form a Bloch 
wall (or soliton) . Such Bloch walls may have various ori- 
gins in realistic samples. They can simply be enforced 
by keeping the spins at both sample ends antiparallel to 
each other. For certain sample geometries, their existence 
can be favored by long range m agne tostatic interactions 
which have not been built into (2.5). Finally, in strictly 
one-dimensional chains, solitons with width of one lattice 
constant rathjet-than spin waves can form the elementary 
excitations.EjH 

A static Bloch wall connects the anisotropy minima 
= ±^ within the easy plane = f and thus satisfies 
the Eulcr-Lagrange equations 



Jdi 



Ky sin ( 



<cos(h = 0. 



(3.1) 



With the additional condition dx 
immediately integrated once 

-^(9,0)2 -cos2<^ = 0. 

Ky 

This equation exhibits the symmetries 
(j) ^ (f) + TT, which reflect the fact that the energy ( 2.1 ) is 
invariant under rotations by tt around each axis in spin 
space. Cons eque ntly, there are four different Bloch wall 
solutions of (Ol) (see, e.g., Ref. E^) 



(3.2) 



and 



(j>Qc{x) = -QC- 2 arctane'^^/^ 



f, (3.3) 



of width 5 = yJjjKy. In order to distinguish the four 
different soliton configurations we have introduced the 
"charge" Q — ^ Jdx dx{sm(j)), and the "chirality" 



dx dx(i 



(3.4) 



of a spin configuration. For the Bloch walls (3^) we have 
Q. C — ±1, and all four walls have the same energy 

/oo 
dx {dx^Qcf = ANa^/JK^. (3.5) 
-oo 

The definition of the chirality C simply tells us whether 
the angle (j) increases or decreases as we proceed in pos- 
itive x-direction along the sample. The definition of 
charge Q is motivated by the response of the Bloch wall 
to an applied magnetic field: For an external field along 
the positive y-axis, a Bloch wall of positive charge moves 
along the positive a;-axis while a negatively charged wall 
moves in the opposite direction. [We recall that the spin 
is antiparallel to the magnetization.] Within the present 
description, the spin is allowed to point into an arbitrary 
direction on the unit sphere . In this case, only the 
charge is a topological invariant, i.e., for infinite sample 
length field configurations of opposite charge cannot be 
deformed into each other without overcoming an infinite 
energy barrier. Solitons of different chirality (but same 
charge) can be deformed into each other via a "Neel-wall" 
configuration where the spin at the wall center points 
along the hard-axis. It is only in the ATy-limit of large 
hard-axis anisotropy where the configurations space of 
the spins becomes a circle and the chirality also becomes 
a topological invariant. It will be this limit which shall 
be considered in the next section, but we shall return to 
the gene ral case when we discuss chirality tunneling in 
Sec. IvTl, 



IV. RELATION TO THE SINE-GORDON MODEL 

In some materials such as elongated YIG films (cf. Fig. 
^ or in garnet crystals with misoriented anisotropy axes, 
the hard-axis aniotropy is much larger than the easy-axis 



4 



anisotropy, typically by a factor of ten or more. As a con- 
sequence, deviations away from the easy-plane become 
energetically costly and the magnetization will be con- 
fined to the easy-plane and the system can effectively be 
described in the easy-plane variable cj) only. 

In the limit Kz 3> A'j,, deviations awav-^rom the easy- 
plane are suppressed and we can expands 



(4.1) 



where |t?| <C 1. Inserting (4J) into the action (|2.4|) we 
obtain up to second order in 'd, 

Se = Na Jdxdri^ i —drcf) + J{dx(t)f + KyCos^ 



i-d dr(l) + ^ L d 
a 



(4.2) 



where L = -Jdl - J{d^(j)f + Ky sin^ (j^ + K,,. If the fluc- 
tuations in both ■& and (p have wavelength A larger than 
the domain wall width, A > (5, the hard-axis anisotropy 
becomes dominant and C = + 0{ Ky/K z)- With 
this approximation, we insert (4.1) into ( |2.2D and using 
VcosO ~ V-d we can perform the Gaussian integrations. 
The transition amplitude can then be expressed as a path 
integral over the azimuthal angle (p alone 



{{n,}\e-^^\{n,}) 



V(j> e 



-Ssg[ 



(4.3) 



with the boundary conditions (f){x,0) — (j)a{x), (j){x, P) 
(j)b{x). The action has the following form 



a 



1 



SSG = Na ldxdT{ i-drP+J[^{dr^f + (3,0)2] 



+ Ky COS 



(4.4) 



where we have introduced the asymptotic spin wave ve- 
locity. 



(4.5) 



We thus have arrived at the important result that for 
large hard-axis anisotropy, the dynamics of a mesoscopic 
ferromagnet is described by the sine-Gordon (SG) action 
plus a topological term iSNA J^Jdrcj). While the re- 
duction to— the sine-Gordon model has been known for 
some tima^J'El, the topological term has not been identi- 
fied before. This term is of central importance for the 
quantization of the spin system as we shall see below. It 
is this term that is responsible for observable effects such 
as band halving and chirality correlation. 

We can now explicitly verify the consistency of our ap- 
proach. SsG has the same long wavelength excita tion s 
as the full magnetic model described by the action (2.3). 



In the latter model, deviations from the uniform state 
(f) = ^,6 = ^ along the easy-axis have the spin wave spec- 
trum ujk = 2{a/S){[JP + Kz][Jk^ + Ky]^/^, with 
k the spin- wave wavevector£3 For ^ Ky and k < 
v/iWJ this reduces to uJk = 2{a/ S){Kz[Jk^ + Ky]^/^ = 
c\Jk^ -|- (5^^ which is identic al t o the spin wave spectrum 
in the sine-Gordon model (4.4) around — 11/2. Sim- 



ilarly, the dynamic soliton solutions of the spin system 
(see, e.g. Rcf. ^) which correspond to moving Bloch 
walls have their counterpart in the SG-model in this 
limit. Even soliton- antisoliton breather solutions of the 
spin system have analogues in SG breather solutions.EZI 
This is surprising since in the spin model breathers ex- 
hibit a precession around the easy-axis and thus do not 
stay close to the easy plane as required in (4.1). The 
connection between the spin model and the SG system 
is therefore more general than the above derivation sug- 
gests. 



A. Spin Tunneling for K-^ ^ Ky 

To illustrate the importance of the topological term 
derived above, we consider the case of a uniform spin 
configuration as, e.g., realized in a nanoscale ferromag- 
netic particle. We shall show that the reduced model 
(4.4) reproduces^ both the spin parity effectli3 and the 
tunneling actionO of the full magnetic model in the XY- 
limit. In addition we shall also evaluate the prefactor of 
the transition amplitude resulting from Gaussian fluctu- 
ations around the instanton path. 

For uniform configurations, <f) = (/'(t), the action (4.4) 
reduces to 



SsG = N dr 



j dT^iSdr(l} + ^{drpf + Kyacoi? (1)^ 



(4.6) 



where N = NaL/o is the total number of spins in the 
sample. Note that </> describes the azimuthal angle of 
a spin and is a compact variable {(f) -\- 2it is identified 
with (/)) . The tunneling amplitude between the anistropy 
minima at = ±7r/2 is then given by 



-f3H\ 



2' 



0(/9) = 
0(0) = - 



e 



-•Ssg[' 



(4.7) 



The dominant contributions to the transition amplitude 
are the extrema of the action which satisfy (JiSsg = or 



Ky sin ( 



I cos (i = . 



(4.8) 



Note that this equation is formally equivalent to ([3.l|). 
Similarly, as a consequence of the symmetry of the ac- 
tion (4.4) under 7r-rotations around the hard-axis, i.e. 
^ 4> -\-TT, the two anisotropy minima can be connected 
by two different paths. These "instanton" and "anti- 
instanton" trajectories are given by 4>±{t) = —tt/2 ± 
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2arctane"^*^'^~'^°-' and describe a transition from = 
—■k/2 to 7r/2 in clockwise (</>+) or anticlockwise direc- 
tion (0-). The transition occurs at tq within a finite 
imaginary time interval characterized by the "instanton 
frequency" w/ = c/ 5 — {2a/ S)yJ KyKz- Inserting (f)± 



into the action (4.6), we recognize that the topological 
term gives rise to a phase which differs in sign for instan- 
tons or antiinstantons 



'5sg[0±] = ±iTTNS + So, 



(4.9) 



with the tunneling exponentEa Sq = 2N S Ky/ Kz- 

The effect of the topological phase may now be seen as 
follows. Adding the contributions of one single inst anton 
and antiinstanton to the action, we obtain for ( [4.7D 



10: 



1\ 
2' 



cx e 

cr=±l 

2cos(7riV5)e"'^° 



Ssg[ 



(4.10) 



Thus the transition amplitude vanishes for half-odd 
integer NS since tunneling paths of opposite wind- 
ing (or "chirality") interfere destructively with each 
other. A calculation within the "dilute instanton gas" 
approximationO reveals that this interference persists 
to all orders in the instanton contributions. Identify- 
ing /i = 2NJa/c^ = NS'^/2Kza, k = NKya, a = -NS, 
d — n we obtain (including the contributions of Gaussian 
fluctuations around the instanton) 



-/3-H| 



with A = 16 



/? ^ oo in Eq. (4.11 



sinh ( cos(iVS'7r) ) . (4.11) 



a (K^Kz) e-^° . Taking the limit 
we conclude that the ground state 



energy of the individual potential wells uJi/2 is split into 
two levels separated by A provided the total spin NS 
of the particle is integer. For NS half-odd integer no 
splitting occurs. Thus for arbitrary spin S the splitting 
AE between the two states of lowest energy is given by 



AE = |cos(7rA^S')|A. 



(4.12) 



Thus, in the uniform limit our theory thus reproduces 
the spin-parity effect of Ref. moreover, in the limit 
Ky ^ Kz, the tunneling exponent agrees with Refs. Igc" 
and A with Ref. |. 



V. BLOCK WALLS AND SPIN WAVES 

In this section we discuss the interaction between a 
Bloch wall and its surrounding spin waves. We consider 
here a sample with Na (or S) sufficiently large such that 
spin waves are just a small perturbation of the Bloch wall. 
For a quantitative description of this interaction we use 
a systematic approach with the ratio of wall velocity to 
spin wave velocity X/c as a small parameter. This is 



justified since typically X <^ c lO^cm/s. We con- 
struct then an ab initio theory for the soliton dissipation 
by integrating out the spin waves. Finally, by deriving 
the spectral function of the damping kernel we can make 
contact with the phcnpmenological Caldeira-Leggett for- 
malism of dissipation.t3 A brief account of the following 
results has been given in Refs. [2l|j2^ . 

We consider elongated samples (cf. Fig. ||) of suffi- 
ciently smalLpipss sectional area such that the transversa 
spin wavesEjEj around the Bloch wall are frozen out.c£l 
This condition typically requires N^ to be less than 10"^, 
and thus can easily be reconciled with above condition 
that Na ^ 1. Motivated by materials such as YIG which 
are favorable for MQP we consider the limit of large hard- 
axis anisotropy. This allows us to build upon the results 
of the last section and we can treat the interaction be- 
tween Bloch wall and spin waves within the sine-Gordon 
model.E^ 

For notational simplicity we restrict ourselves for the 
moment to one of the Bloch walls ( |3.3[ ), (j)o{x) — 
4'Q=i,c=i{x)- First we recall that (j)Q{x — X) is, for arbi- 
trary X, a static solution of SSsg ~ 0. We now consider 
field configurations describing a Bloch wall at a position 
X surrounded by arbitrary spin waves ip 



(x, t) = 0o(a; - A) + Lp(x - A, r), 



(5.1) 



and elevate A(r) to a dynamical variable. However, Eq. 
( ^.1[ ) contains now a redundant description of a rigid 
translation of the soliton: A translation is either de- 
scribed by X or|-bjy|-tJie, "zero- mode" (Goldstone mode) 
(y9o(a;,r) cx 0Q(a;)l^ c3'Ej. To avoid double counting, we 
thus have to impose the constraint that the spin wave 
modes be orthogonal to the zero mode 



dx (/)g(x)(p(x, r) — 0, 



(5.2) 



for all imaginary times r. We incorporate this constraint 
into the Mtiii integral by means of the Faddeev-Popov 
techniqueoS which we now briefly sketch. It is based 
on the identity 



VX5(Q[X])det 



5Q_ 

SX 



1, 



with the judiciously0'0 chosen functional 
Q[X] = / dxq)'^{x~X)(t){x,T). 



(5.3) 



(5.4) 



Inserting (5.1) into ( |S.4[ ) we recognize that the (5- function 
enforces the constraint ( ^.2[) as desired. 

For configurations which contain o ne s oliton, we thus 
can rewrite the transition amplitude ( [4.3| ) as follows 



{{nt}\e-^^\{na}) = JvXVcj, 5{Q) det|| e-^-'^I 



(5.5) 
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where the action is given by (4.4). We now perform a 
systematic expansion up to 2"^^ order in both spinwaves 
(p and X/c. [Note that X/c < 1.5 • 10"^ foj. yiG as 



discussed in Sec. VII and Lp o c 1/\/ Na as we shah see 
befow]. After insertion of (5.1) into (5.5) and expansion 



to second order in the spin waves ip and second order 
in the Bfoch wall velocity X/c, the transition amplitude 
takes the form 

{{nh}\e~^^\{na}) ^ fvX e-^- W F[Xl (5.6) 



where 



Sx = I dr|-mX+yX2| 



is the action of a free Bloch wall, and where 
F[X] = Jv^ ,5(^/0;,^^ det If e-^^{^-[e+'Cl^+'^-^> 



(5.7) 



(5.8) 



describes the interaction between the Bloch wall and the 
spin waves. Here we have introduced the scalar product 
a ■ b = J dx dr a*b and the integral in the (5- function is 
understood as an integral over x. 

We now discuss the various terms that have been in- 
). The first term in this action has 



troduced in (|5.6|)-(|5 



the form of a gauge potential 

T:SNA/a. 



(5.9) 



It originates from the topological term in (4^) and from 
the relation / dx dr (t){){x — X) = —ttX since each soliton 
flips the spins by tt = / dx^Q . 

The second term in ( [3.7] ) is the kinetic energy of the 
Bloch wall and the mass is given by 



M 



En 




(5.10) 



This value coincides with the Doring mass0 M = 
7Vyia^/(27r7^(5) with 7 = gfis/^ if the hard-axis 
anisotropy is of purely demagnetizing origin, = 
27rMQa^ with Mq — g^sS/a^. We thus have given a 
microscopic derivation of the Doring wall mass. 

In (5.7) we have dropped a term PEq with Eq the Bloch 
wall energy (3.5) since the Bloch wall already exists in 
the sample and is not created thermally. The thermal 
creation of Bloch wall pairs in the absence of an ex- 
ternal field is negligibly small for temperatures in the 
Kelvin range even for samples as small as 50A x 50A. 
Only at higher temperatures and in the presence of ex- 
ternal fields -thermal creation of Bloch wall pairs becomes 
appreciable .E3 

The functional F ( |5.§| ) describes the couplijaff between 
spin waves and the Bloch wall. The operatoio 



Ky[l 



2sech2(^)] 




(5.11) 



with K = J /(? describes the spin wave spectrum around 
a static Bloch wall. The remaining operators are respon- 
sible for the dynamic coupling between spin waves and 
domain wall. 



K, = 2KXd^dr~ K.X'dt J 



-2KX'^(t>Q. 



(5.12) 



Due to the constraint (5.2) the exponential in F, Eq. 
( ^.SD , does not contain a term linear in the velority X 
and in the spin waves ip. It has been pointed outcHlEil that 
this is an important difference to the standard Caldeira- 
Leggett modeLcJ However, despite this non-linear cou- 
pling we shall see shortly that at low temperatures the 
dissipation due to spin waves can — if this should be 
desirable — perfectly well be modeled by a Caldeira- 
Leggett approach (such cases have actually been dis- 
cussed in App. I of Ref. ^5|), although the precise form of 
the relevant spectral function can only be obtained from 
a microscopic calculation as presented here. 

After the evaluation of ( |5.5| ) w hich is sketched in Ap- 
pendix ^, and collecting Eqs. (5.7), (B4), ( p5|), and 
(B17) we can express the transition amplitude (|5.6|) as 



{{n,}\e-^^\{na]) = VXe 



-Sett[X] 



with the effective action for the soliton position. 



(5.13) 



/3 



drl —iaX 



M, 



off 



X^} 



)?■ (5.14) 



7: [ dr [ daK{T-a)[XiT)~X{a 
2 Jo Jo 



The damping kernel has been evaluated for arbitrary 
temperatures in Appendix ^ Here we restrict ourselves 
t o low temperatures, /? — > 00, where the damping kernel 
(B19) takes the following form 



-2wfc|r| 



(5.15) 



This cai^alao be cast into standard Caldeira-Leggett 
notatiorMEa 



1 r 

K[t) = - / 

Jo 



duj J{uj)Di^{t). 



(5.16) 



Here D^{t) ^ e^"^!^! is the T-^ hmit of ( |B1^ ) and the 
spectral function is given byE2 



^0 ' 



(5.17) 



which vanishes for w < cjq = 2c/5 — {Aa/ S)y/KyKz, the 
anisotropy gap of the spin waves. [For material values as 
in YIG, this gap corresponds to a temperature of Tg = 
0.2K. Other materials have in general larger anisotropics 
and thus higher Tg.] In derivin g (^.171) from (5.15) we 
have used the renormalization (piO[). It is only after 
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this renormalization that the memory kernel K becomes 
positive definite (as is needed for convergence). 

If now the dynamics of X is slow compared to the time 
variation of the damping kernel, i.e. if the instanton fre- 
quency 07/ (to be evaluated below) is much smaller than 
luq — 2c/S, and if the temperature is small such that /? ^ 
LuJ^ then we may expand X{t) — X{a) « (r — a)X{(7), 
and the damping kernel reduces to a pure mass renor- 
malization. Note that this mass renormalization is a 
0{{Na)'^) correction of the wall mass M oc Na- Since 
all these conditions will be satisfied for the tunneling sit- 
uations considered below, we see that we end up with a 
deceptively simple effective description of the B loch wall 
dynamics, given by the first two terms in (5.14). 



VI. INTERFERENCE EFFECTS DUE TO THE 
BERRY PHASE 

In the last section we derived an effective action for 
the dynamics of the Bloch wall position. We showed that 
damping due to spin waves leads to a gap in the spectral 
function and thus leads to a mere renormalization of the 
wall mass at low temperatures. More importantly, we 
have identified a topological term in the action which 
has its origin in the Berry phase term (2.4) of the original 
spin action. 

Here we shall generalize these considerations to soli- 
tons 4>QC of arbitrary chirality C and charge Q which 
are all energetically degenerate. We show that the in- 
terference between states of different chirality gives rise 
to remarkable effects such as the halving of the Brillouin 
zone and the alternation of chirality in reciprocal space 
for half-integral spin. 

These effects originate in t he f act that for arbitrary 
solitons, the topological term (2.4) 



i— I dx(pQc{x — X) = —iaCX, 



(6.1) 



depends on the soliton chirality C (but not on its charge 
Q) (a = ttSNa/o,). This chirality dependence can intu- 
itively be understood as follows. As the soliton sweeps 
across a given spin, the spin is rotated by an angle ±7r 
(after the wall is sufficiently far away), the rotation sense 
being uniquely determined by the chirality and the di- 
rection of motion of the Bloch wall. 

At low temperatures a nd fo r solitons of arbitrary chi- 
rality the effective action ( 5.14 ) of a soliton in an external 
potential V{X) thus takes the following form 

S[X,C]^ J dT\^~iaCX + ^X^ + V{X)y (6.2) 



Here we have used the mass ( 5.10 ) rather than the 
dressed mass M^s since the mass renormalizations at low 
temperatures are small 0{{Na)'^) and the value ( S.lOp 
thus represents a good approximation for the experimen- 
tally observed wall mass. In addition we have introduced 



a periodic potential V{X) of period d, and we make the 
natural assumption that d is some integer multiple of the 
lattice constant a. For definiteness we assume 



F(X) = yo[l-cos(^^)] 



(6.3) 



which has-ipnplitude 2Vb.0 Such a potential can have 
its originoO in the discrete nature of the crystal lattice 
itself, or for Bloch walls it can arise from a magnetic 
superlattice of lay ers with different anisotropics. 
The action (|6.2|) corresponds to the Hamiltonian 



n 



aa,f + V{X), 



(6.4) 



where p — —id/dX is the Bloch wall momentum, az 
is the Pauli-matrix of the "pseudospin" characterizing 
the chirality C = ± of the Bloch wall. Obviously, this 
Hamiltonian conserves the chirality. For mathematical 
convenience, we choose periodic boundary conditions in 
the following. However, all our results are finite in the 
thermodynamic limit and none of our cpmclusions depend 
on this choice of boundary conditions £3 

From both ( |6.2[ ) and (6.4) it is evident that the topo- 
logical phase plays the role of a gauge potential whose 
effect on the wall dynamics shall be discussed next. We 
note that such spin-dependent gauge potentials-are not 
uncommon in problems involving Berry phases.tJ 



A. An illustrative example 

Before giving a rigorous discussion of the dispersion 
relation, we give an argument to illustrate the interplay 
between the topological phase and soliton propagation. 

Consider the transition amplitude for the propagation 
of the Bloch wall between nearest neighbors which is 
given by 



c -^0 

= 2cos(ad) / VXe-^°, (6.5) 
Jo 

where 5o[X] = J dT{{M/2)X^ + ViX)}. For half-integer 
a = NASd/a, we thus arrive at a most important con- 
clusion: Nearest neighbor hopping of the soliton is sup- 
pressed if both chiralities contribute equally to the tran- 
sition amplitude. However, if the soliton is in a state of 
definite chirality, only one path contributes to the tran- 
sition amplitude (6.5) and nearest neighbor hopping is 
allowed. No such interference occurs for integer a. Note 
that this interference effect is entirely due to the topolog- 
ical term in (6^) which in turn is a conseq uen ce of the 
topological term in the sine-Gordon action (4.4). 

We now investigate how this interference affects the 
dispersion of solitons and in this way can become ob- 
servable. 
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B. Dispersion in the nearly free limit 

In this section we discuss the dynamics of a sohton in 
an arbitrarily weak periodic potential V{X). The Hamil- 
tonian is given by (6.4) with Vq ^ 0. Despite being sim- 
ple this case already captures most of the characteristic 
features of the tight binding limit which shall be dis- 
cussed below. For simplicity we assume that the period 
of the potential is given hy d — a. 

U sing periodic boundary conditions, the eigenstates of 
(3.4) are simply plane waves e'*^"^, with k — 2TTn/L, 
L — Nd, and the spectrum consists of two parabolas 
(corresponding to the two soliton-chiralities) 



E{k,C = ±l) = ^{kTaf, 



(6.6) 



periodically extended by the reciprocal lattice vector 
G — 27r/a. [Note that the requirement of gauge invari- 
ance alone produces such a periodic extension even in 
com plete absence of a perio dic potential. The gauges 
(A3) lead to Hamiltonians (6^) with a (2n + l)a 
with n = 0,±1,.... The gauge in vari ant dispersion is 
therefore the periodic extension of (|6.6| ) by a vector 2a.] 
The Berry phase thus leads to remarkable spin parity 
effects in the dispersion: For half-odd integer spin NaS, 
we have a — G/A (mod G) and the parabolas are sep- 
arated by half the reciprocal lattice vector, G/2. Thus, 
the Brillouin zone is halved and two subsequent parabo- 
las belong to opposite chiralities as illustrat ed in FigJ|. 



The observability of this is discussed in Sec. VI D 



For integer spin, however, the dispersion is analogous 
to that of a free particle of mass M and the Berry phase 
is inoperative since it merely shifts the dispersion by a 
reciprocal lattice vector. 

Note that since a oc 1/a is independent of the sample 
length, the result is unchang ed i f we pass to the ther- 
modynamic limit. Therefore (6. 
boundary conditions.L3 



is independent of the 



C. Tight binding limit 

We now turn to a discussion of the system in the tight 
binding limit where V{X) is no longer small. In the ab- 
sence of tunneling there exists a large number of degener- 
ate ground states corresponding to the soliton trapped at 
one particular pinning site. If the pinning potential is not 
too large, the soliton can tunnel between the sites, and 
these ground states split into the (lowest) band E{k, C) 
with 



n\k,G) = E{k,C)\k,G). 



•7) 



Since Ti in (|6.4[) is invariant under translations by the po- 
tential period d and conserves the chirality pseudospin, 
the eigenstates are products of Bloch states and chiral- 
ity eigenstates C) = |fc) |G) where Tjfc) = e''"^\k), 



a^\C) = C\C) with fc = 27rn/iVd, n = 0, 1, iV - 1, and 
L = Nd. T is the translation operator {X\T ^ {X + d\. 

For the evaluation of the band structure in the tight- 
binding limit we now develop a formalism which allows 
us to keep carefully track of the topological phases within 
the instanton approach. To-.this end we start from the 
modified partition functionO 



Z,=tr{T'e-'^«}, 



(6.8) 



where tr{-} = Efc,c(^'C'l ' l^'^*)- We use (|j) rather 
than the usual partition function Z = tr{e~'^'^} for the 
following reason. As we have seen in the previous section, 
the Berry phase gives rise to a shift of the dispersion with 
respect to k. However, the partition function Z is insen- 
sitive to such shifts (at least if there is no perturbation 
which mixes the chirality states) and thus represents an 
insufficient tool for the evaluation of the band structure. 

From ( ^.^ ) we can easily extract the dispersion by tak- 
ing the Fourier transform 



AT-l 

(=0 



-(3E{k,C) 



(6.9) 



c=± 



where we used (6.7) and the definition of T. In addition, 
we have restricted ourselves to the lowest band since we 
are interested in the low temperature limit. Note also 
that the Ihs of ( |6.9| ) does not contain higher winding 
number contributions since we are not interested in fi- 
nite size effects arising from the sample topology. We 
now evaluate Zi and get in a first step, 

Zi=y2 [ dX {X + ld,G\e-^'^\X,G) 

N-l f.md+d/2 

= E E / '^^ I G|fc) |2e-'3^^('='t^)+^'=''^ , (6.10) 

C,k m=0 '^/2 

where we used periodic boundary conditions and inserted 
a complete set of Bloch states. Next, in the tight binding 
limit the main contributions to the integral are coming 
from the vicinity of the potential minima, X ~ md, and 
the Bloch functions can be replaced by their harmonic 
approximations, i.e. \{k\md)\'^ « |V'h(0)P/A^ = a^/N. 
Here, V'h is the ground state in the harmonic approxima- 
tion of the potential well and oq — ^ Muj / vr its normal- 
ization (squared). Thus, we find 

N-l 

ZiK — ^^{md + Id, G\e-'^^ \md, G) . (6.11) 



C m=0 



Using a path integral representation for ( ^.ll }_and em- 
ploying the periodicity of Ti. we obtain with (6.9) 



E' 



-f3E{k,C) 



ao 



E' 

C,l=0 



ikld 



X{p)=ld 



X(0)=Q 



VX e-^^^'^l (6.12) 
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The path integral on the rhs of (6.12) is dominated by 
instantons between the potential minima. These in- 
stantons obey the Euler-Lagrange equation 5S/5X — 
—MX + V'{X) — 0. For instance, a transition from 
X = {)ioX — d{X — ~d) is mediated by the (anti-) 
instanton 



X^ 



2d 

± — arctan e' 

TT 



(6.13) 



centered at the arbitary imaginary time tq. The instan- 
ton frequency uj — ^a/Vq/M equals the harmonic os- 
cillation frequency in the potential well. The instanton 
action is given by 



(6.14) 



where S y, — ^ d\/MVo — %[Vo/lo). The miusual second 
term in ( 6.14 ) is purely imaginary and is a dire ct conse- 



quence of the gauge potential in (5.14) or (S.4) and dis- 
tinguishes between instantons and antiinstantons. Note 
that this term does not break time reversal invariance 
as the partition function contains contributions of both 
chirality states C — ±1. 



The path integral in ( 6.12 ) can be expressed as the sum 
over all distinct sequences of instantons and n_ — 
n+ — I antiinstantons which connect the initial state X =^ 
with the final state X{[3) ^Id. Within this "dilute 
instanton gas approximation" t3 we obtain 



X(P)=ld tx) 

VX e-5[^^^l = floe-- X 

-^(0)=0 n+=0,n-=0 



{JK[3e 



-S+\nj, 



{JKI3e 



-S- \n- 



(6.15) 



where J — Sq/2ttM ^ K — 2a;vM arise from the in- 
tegration over the zero modes and the Gaussian fluctua- 
tions around an instanton, respectively. Inserting ( 6.15| ) 
into ( |6.12 ), using ( |6.14 ) and performing the sums we ob- 
tain 



E 



G 



where 



'2JKe~^" cos 



i(^{k + aC)d 



The ground state as a function of k is given by 
1 



lim 



c 



(6.16) 



(6.17) 



(6.18) 



Similarly to the nearly free limit discussed in the previous 
subsection, this dispersion is fundamentally diffepent for 
a = ad/n — NASd/a in teger or half-odd integer .Eil 

Inserting ( 6.16| ) into ( 6.18 ) we obtain for integer a the 
following dispersion 



E{k) = -— cos(fcd), 



(6.19) 



which is of standard tight-binding type. In (6.19) we 
dropped the constant uj/2. The bandwidth is given by 



(6.20) 




In contrast, for a half-integer we obtain 



A 



sin(fcc?)|. 



(6.21) 



In (6.19) and ( |6.21 ) we have suppressed sign changes 
which correspond to a global shift of k by ir/d. Such a 
global sign cannot be measured since the absolute value 
of k is experiment ally n ot detectable. 

The dispersion ( 6.21 ) has now cusps and the bandwidth 
and the Brillouin zone a re ha lvecE^, as s hown in FigJJ. 
Moreover, we draw from (|6.16| ) and ( |6.18 ) the important 
conclusion that states whose wavevector differs by tt have 
opposite chirality, cf. Fig.^ 

Note that this period halving in reciprocal space is a 
consequence of t he fa ct that Zi = /p'' 'DXe~'^^^''~'^ = 
for / odd (cf. (6.5)). However, one must not conclude 



from this fact that nearest neighbor hopping is always 
supp resse d: At fixed k ^ 0, ±-|, the ground state condi- 
tion (3.1S) selects a branch of the dispersion with definite 
chirality, a dispersion that arises from nearest neighbor 
hopping. Only at the cusps in ( 6.21 ) hopping is sup- 
pressed. 

This band-halving can also be understood in a more 
intuitive way: For half- integer a, a soliton acquires a 
Berry phase —iC for forward {iC for backward) hopping. 
In the ground state this phase gets compensated by the 
Bloch phase thus creating two band minima at kd = 
±7r/2 which have opposite chirality. 

Finally we give a more explicit formal argument for 
the chirality correlation. We find the explicit form of the 
eigenv alue E{k, C) by repeating the steps leading to Eq. 
il), 



, but instead of Z/ we use = tr{|C)(C|r'e-'3«}, 

which projects onto a state of definite chirality C. T hus, 
we find that E{k, C) is given by e(fc, C) in Eg. ( |6.17| ). By 
comparing the ground state energy ( |6.21 ) with E{k,C) 
we see that fc-intervals with positive (negative) sin kd be- 
long to negative (positive) chirality C. This result is de- 
rived in the north pole parametrization. If, instead, we 
use the south pole parametrization, then the gauge po- 
tential in (S.2) changes sign and again we find that the 
chirality alternates, but now with the opposite assign- 
ment between chirality and given fc-interval. The physi- 
cal consequence — alternating chirality with changing k — 
is the same in the two gauges, since, again, the absolute 
k value cannot be observed. 
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D. Discussion and analogies to other physical 
systems 

In the last two sections we have seen that the dispersion 
is strongly affected by the parity of cr = NASd/a. For 
(T integer, the dispersion equals that of a particle in a 
periodic potential while for a half integer a halving of the 
Brillouin zone occurs with alternating chiralities. In the 
latter case the dispersion consists of mutually intersecting 
parabolas or tight binding bands. How can we observe 
such a dispersion? 

Let us for definiteness focus on the nearly free limit 
with a dispersion as shown in Fig. |^. Suppose the chiral- 
ity has been measured to be C = 1 and the system is in its 
ground state, i.e. in the minimum of a C = 1 parabola. 
If we now drive the system out of its energy minimum, 
e.g., by applying an external field along the easy-axis (see 
below), the Bloch wall will follow the C — 1 parabola. 
The Bloch wall will remain on this parabola even be- 
yond the crossing point, provided that the chirality C is 
a conserved quantity. In this sense, the two parabolas for 
C = 1 and C = — 1 behave like two "Riemann sheets" of 
the energy which are completely disconnected, and their 
intersection has no observable consequences if there is no 
mixing i.e. tunneling between the chiralities of the Bloch 
wall. 

Nevertheless, the dispersion of Fig. ^ (thin line) is a 
precursor of a striking physical effect: As soon as there is 
tunneling between the chiralities, the different "sheets" 
get connected and for half integer spins a gap develops 
at the crossing points of parabolas belonging to C = 
±1. At the same time the halving of the Brillouin zone 
becomes observable. Formally this can be described as 
follows. In the prese nce o f tunneling between the two wall 
chiralities (see Sec. VII) the Hamiltonian ( |6^ ) acquires 
an additional term ea^ 



n 



2M 



(6.22) 



such that the chirality C (i.e. Uz) is no longer a conserved 
quantity. We are interested in the limit of small chirality 
tunneling and therefore e will be much smaller than the 
bandwidth A (estimates for e will be given in Sec. VII). 
For cr half-odd integer, the degeneracy at the points fc„ = 
WK / a is lifted and the dispersion splits into two bands 
which for |/c — fc„| <C 7r/a are given by 



E±(k) = — 



knf + a' ± v/4a2(fc-fc„)2-f e 



22 



(6.23) 



with e — 2Me and where, for simplicity, we have stated 
the result in the nearly free limit. In this and the tight 
binding limit the two bands are separated by a gap 2e at 
k = kn, as shown in Figs. ^, |. 

Solving for the corresponding eigenstates we recognize 
that the chirality continuously switches from C = ±1 



to C = ^1 as we pass from one band minimum to an 
adjacent one. 

W e t hus have established that the spectrum given in 
( |6.6|) , ( 3.21 ) is reached in the limit e — > 0. Note that 
the experimental observation of the gap depends on the 
probability of Zcncr intcrband transitions and thus on 
the time scale at which the band structure is. probed. In 
the nearly free limit, the Zener probabilityO can be ex- 
pressed as0'EI P cx exp{— 11^}, where T ~ 2tt/ujb is 
the time to cross the Brillouin zone, with lob = Fd/h be- 
ing the Bloch frequency and F = 2g^BSNAH/a the driv- 
ing force due to an external field H (along, say, the easy 
axis, see Sec. VIII). Eq ~ j^{'K/d)'^ is the kinetic energy 
at the zone boundary. Thus, to optimize observability 
we must have A — ^j-p- ^ 1, which is easy to achieve 



since typically T 
choose H 
(see Sec. 



VII) 



2 hEo 

IQ-'^s, giving A - 100 for YIG, if we 
lO-^Oe, d = a, and e - Eq/IO - lOmK/cs 
The alternation of chirality could then 



be observed, for instance, by (optical) dichroism tech- 
niques which would be sensitive to the rotation sense of 
the magnetization within the Bloch wall. 

We emphasize that these results are gauge indepen- 
dent. If, instead, we had started from the south pole 
parametrization of the coher ent st ates, w e wou ld have 
obtained the same dispersion (|6j), (|1|), (|^, except 
for a global shift k k + 2a which is unobservable. 

A dispersion consisting of disjoint parabolas dictated 
by gauge invariance and the formation of gaps due to tun- 
neling is quite a common phenomenon in condensed mat- 
ter physics. Persistept currents in isolated metal ringsO, 
the Josephson effeclllj, and the tunneling of quasiparti- 
cles beiween edge states in the fractional quantum Hall 
regimeES might serve as familiar examples. 

For further illustration let us briefly discuss some rela- 
tions between our spin effect and, say, persistent currents. 
First, in the spin system the dispersion remains unaltered 
in the thermodynamic limit, whereas persistent currents 
are a finite size effect, resulting from the discreteness of 
the energy levels. In addition, we consider a simply con- 
nected sample topology while a persistent current relies 
on the ring geometry of the sample. In the spin system 
it is the 5^-topology of spin space restricted to the easy- 
plane, not the topology of the sample which is responsible 
for the interference effect. 

An electron of mass m confined to a ring of radius g 
which is threaded by the electromagnetic fiux $ is de- 
scribed by the Hamiltonian 



n 



2niQ' 



(6.24) 



where $ is measured in units of the fiux quantum $0 = 
/ic/e, and 9 is the azimuthal angle. The eigenfunc- 
tions are e™* with eigenvalues £"„ — (fi^ /2mQ'^){n — , 
where n = 0,±1,... The ground state energy Eg as 
a function of flux is the ens^lope of the set of energy 
parabolas separated by ^o-l^ Thus, the persistent cur- 
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rent j = —{e/h)dEc/d^ is a sawtooth curve with dis- 
continuities at |<i>| = n/2 where the parabolas intersect. 
Suppose now that $ = and that the system is in its 
ground state with n — 0. If the flux is increased adiabat- 
icahy, the system will stay on the n — parabola even 
for $ > 1/2 since the angular momentum is a conserved 
quantity. Thus the electron will not see the other parabo- 
las and the spectrum consists of disconnected "sheets" of 
parabolas. This b eha vior is analogous to that for the 
soliton dispersion (6.6) for half-integer spin. 

However, if angular momentum is no longer conserved, 
e.g., due to the presence of a scattering potential, the 
parabolas will be connected and a gap develops at their 
crossing points. The scatteri ng po tential thus plays a 
role similar to the eax-term in (|6.22|) caused by tunneling 
between the chiralities. 

The mere existence of interference effects in a metal 
ring can also be derived from the following argument. 
Assume that $ = 1/2 and let us imagine having pre- 
pared two wave packets of opposite angular momentum, 
but otherwise identical. If we let these wave packets dy- 
namically evolve until they have traveled half the circum- 
ference, one clockwise the other anticlockwise, they will 
have picked up Aharonov-Bohm phases of opposite sign 
such that (for $ = 1/2) destructive interference occurs 
leading to a vanishing transition amplitude between ini- 
tial and final state. Th is behavior is similar to the spin 
case described in (6^), where the clockwise and anti- 
clockwise traveling wave packets correspond to the two 
chirality states of the soliton (note again that the real 
space topology of the ferromagnet is irrelevant). 



VII. EXPERIMENTAL IMPLICATIONS 

In this section we give numerical estimates for the ef- 
fects discussed in the previous sections. For definite- 
ness we concentrate on material parameters for YIG. 
ExchangeEa and anisotropyEHl are given by J = 1.65 • 
10~^^erg • cm, and Ky = 9.61 • 10~^^erg/cm, where a cell 
with lattice constant a = 6.2A contaiBS,one S = 5/2 spin 
implying a saturation magnetizatioriEa of Mq = 1940e 
(i.e. = = 9.1 • 10-i°erg/cm), wall width 

6 = ^/J/Wy = 414A, and spin wave velocity, Eq. ( [4.5[ ), 
c = 6 • lO^cm/s. The pinning potential strength can 
be related to an experimentally observed coercivity by 
addingpj a Zeeman term — 2^Afo^cxtA to the pinning 



potential V{X), with A = Nao? the cross sectional 
area of the sample. Defining the coercivity He as the 
field at which the barrier height vanishes, we obtain 
Vq/A— HcMod/n. Note that the coercivity is propor- 
tional to the slop e Vp/ d of the potential. Looking at the 
WKB exponent (|T4|), Sq = ^d^/MVo, we see that a 
low coercivity does not necessarily imply a high tunnel- 
ing probability. The crucial condition is a small potential 
width d. |— I 

We now assume a coercivityc3 of He ~ 20e and d — 3a. 



Note that the wall extends over 22 pinning si tes. T he 
instanton frequency then becomes ui = ^y^Vo/M = 

1.4 • 10^°s~^, and \X/c\ = ujid/TTc < 1.5 • 10"^ For a 
sample with cross section A = lO^A we have Na = 260, 
and the wall contains NaS/q = 2 • 10^ spins. The pin- 

SSOmKfcs, 
which is of 



ning potential heigh t tak es the value 2Vb 



and the bandwidth ( |6^ ) is A/h ~ 10 



the order of the measured resonance frequency in Ref. 
[Tsl . The Doring mass, Eq. ( 5.10D , corresponding to this 
cross sectional area A takes the value M = 1.24-lO^^^g — 
1.36-lO^mo, where mo is the electron mass. The crossover 
temperature between quantum tunneling and thermally 
activated behavior is Tc = 2Voh/SokB — tuv/AkB — 
28mK, since Ticu/kB = llOmK for d ~ 3a. Note that 
the bandwidth is extremely sensitive to the details of the 
pinning potential. For instance, ii d = a (lattice pin- 
ning) but all other parameters chosen as above, we ob- 
tain u! = 2.5 • 10^°s~^ (corresponding to 190mK), and 
A/h ~ 1.2 • 10i"s-i (since So/h = 2.3), or A ~ 0.8 
times the pinning potential height 2Vo = llOmKfcs, 
while Tc ~ 48mK. 

We emphasize that these numbers are rather material 
dependent. For instance, in an orthoferrite, a canted an- 
tiferroinagnet, the effective wall mass is by a factor 10'^ 
smalle iO t han the value obtained from the Doring wall 
mass (5.10). Thus tunneling could also occur at much 
larger potential heights and higher crossover tempera- 
tures. |— I 

Next, we briefly address the issue of impuritiesEil, a 
more detailed account will be given elsewhere£j The 
analysis so far was based on the fact that the mag- 
netic field is constant throughout the sample. A sin- 
gle impurity (or similarly an inhomogeneous field) can 
be incorporated into the energy ( ^.5[) by adding a term 
KaS{x—XQ) sin^ cj) where k is of the order of the anisotropy 
constant Ky. Although the impurity is pointlike, it leads 
to an extended potential U{X) = —na sech^((A — xo)/5) 
of width 5 for the wall center. Thus even when na is 
of the order of the strength 2Vo of the periodic poten- 
tial, the impurity potential only leads to a small varia- 
tion {d/5)Ka between pinning sites separated by d <C 5. 
This holds also for a random impurity distribution even 
in the unrealistic case (for YIG) of high disorder with 
one impurity per transverse layer. Under the action of 
an external field Hy along the easy axis, which can be 
much smaller than the coercivity i/c, a-U wells created by 
the impurities can be rendered unstable such that they 
no longer trap the wall. Localization of the wall is then 
determined by quantum intereference effects only which 
we can characterize by the Anderson localization length. 
This length, however, is sufficientiy, large and explicitly 
given by aN\{A/2Vof « 5 • lO^aS 

We note that tunneling in periodic pinning poten- 
tials allows much higher crossover temperatures than 
tunneling out of a single isolated (metastable) poten- 
tial. Indeed, in the presence of an external field 
along the easy-axis the total energy is U{X) = 
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-Vosech^X/S-2AMoH^^tX where Vq = ^ASMoH^. 
The cross-over temperature and the WKB-exponent are 
then given by = 2^/* ^/^fTT^TVlJ^ei/ 



2^/^^nNs^yHJ^TMo€^/'^, where e = l~Hext/H. 



So = 
For ex- 
ample, for a YIG-sample of 50 A x 200 A with He = lOOe, 
this leads to cross-over temperatures in the milli-Kelvin 
range 0.5 mK < Tc < 1.4 mK while the WKB exponent 
changes in the interval 0.2 < So/h < 31.1. 

We now turn to a discussion of quantum tunneling be- 
tween the two chirality states of a soliton. We shall ob- 
tain expl icit e stimates for the level splitting e introduced 
in Eq. ( 6.22| ). In addition, we shall see that chirality 
tunneling provides a novel scenario for mesoscopic quan- 
tum coherence with one important advantage that both 
barrier height and bias of the double well can be tuned 
independently by external fields. 

Chirality tunneling involves rotation of the spins out 
of the easy-plane and thus cannot be described within 
the XK-approximation which we have used so far. To 
treat this case we must go back to the full action ( |2.3| ) 
and deal with both polar angels </> and 0. The gener- 
alization of the wall dynamics to this situation, in par- 
ticular, the reduction to the collective coordinate and 
the dissipation due to spin, waves, is necessarily more 
involved but still feasible .E3 However, since this gener- 
alization is somewhat outside the scope of the present 
work we shall only quote the essential rpaults here and 
give the details in a forthcoming paper.a For definite- 
ness we concentrate now on ferromagnets where the easy- 
axis anisotropy exceeds the one along the hard-axia, i.e. 
Ky ^ Kz] typical examples are bubble materials O To 
take advantage of the resulting approximate symmetry 
around the z-axis, we represent the-jpagnetization field 
as $1 = (sin0sin(/(), cos6',sin0coS(/>).cZl The Bloch wall is 
then described by a rotation of the spins in the a;?/-plane 
about the angle 0, and the chirality switching by a rota- 
tion in the a;z-plane about the angle (p — ±7r. In addition, 
we allow for an external magnetic field along the hard 
axis z with which one can tune the barrier height that 
separates the two wall chiralities. 

Integrating out the 6'-fluctuations around the Bloch 
wall and restricting ourselves to uniform rotations in 
(j) (which is valid£3 if the wall width 5 is less than 
^J J/ Kz),'we obtain an effective Langrangian in 0(r), 



V = Kcos^ -I- 77C0S(/) -I- 77^/4k, 



(7.1) 
(7.2) 



where Mc — Na 



is the effective mass associ- 



ated with the chirality dynamics, and the parameters 
K = 26NaKz and 77 = gfiBSNAT^5Hz/a characterize 
the barrier potential V. Defining the anisotropy field 
by Ha = 'iaKz/ gpLsSi: and noting that the chirality tun- 
nels between the potential minima defined by cos <j)min = 
—Hz/ Ha = u — 1, we obtain for the level splitting e 



where 7 is a numerical constant of order one. The in- 
stanton action Sc and frequency ujc are given by 



V 



(7.4) 
(7.5) 



The crossover temperature becomes — ujc/Sks- Note 
the characteristic power dependence on the external con- 
trol parameter v = 1 — Hz/ Ha with which the chirality 
splitting e can be changed over a large range. In the next 
section we shall also see how a field Hx can be used to 
offset unwanted bias between the potential minima. 

We illustrate these results with some typical numbers. 
Choosing NaS/g - 10^, Ky/Kz ~ 10, 1/ - lO'^, aKz - 
IKfcs, and 5* = 5/2, we find for the chirality splitting 
e w 5mKA:B, while the crossover temperature is Tc ~ 
13mK. The values for the bandwidth A are roughly the 
same as before. This shows that the splitting e can be 
made quite large (on the scale of A) just by tuning the 
external field along the hard axis, while the crossover 
temperature is still reasonably high. Without field, i.e. 
h' — 1, the splitting e is only of non- vanishing value if 
the wall is narrow and/or if Na ~ 1, that means if the 
system is close to being strictly one-dimensional. 



VIII. INFLUENCE OF EXTERNAL FIELDS 

In this section we show that external fields allow us 
to control the gauge potential a. In the presence of ex- 
ternal fields the four degenerate Bloch wall configura- 



tions 



get deformed into new configurations 



0(x), 9{x). For mo ving solitons, (/)(a; 
Berry phase term (2.4) becomes 



5i 



wz 



-iaC I drX, 
'0 



X),e{x-X), the 



(8.1) 



where 



a 



NaS 



+L/2 



L/2 



dx 0'(1 — cos( 



(8.2) 



and the chirality has been defined as C = sgn{/ dx (j)'{l — 
cos 9)} with (j)' = dx4> ■ Note that a is proportional to 
the area on the unit sphere between the north pole and 
the trajectory which is traced out by a given spin upon 
passmg of the Bloch waU (cf. Fig g). Since the Bloch wall 
shape changes in response to an applied external field, a 
will in general differ from the value a — NaStt/u of the 
Bloch wall (pQc, 6* = f- 

An external field is taken into account by adding a Zee- 
man term gi-iB^ ■ X); spin-Hamiltonian Ti. ( ^.l| ). 
Correspondingly, the total energy H (2.5) is changed into 



e = ■i'jLUc\/Sc/2T: e 



(7.3) 
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H = H + NAh ■ dxCl, h = g 



B. 



(8.3) 



For fields along the easy-axis or the hard-axis, the 
static configurations satisfy the Euler-Lagrange equa- 
tions SH = 0, 

Aje' (j)' cos e + 2J(j)" sin + Ky sin sin 2(1) + 
+hx sin (f) — hy cos = 0, 



2J0"+ sin 261 [J(j)'^ - Ky sin^ (/)- 



'hx cos ^ cos ( 



hy cos sin (j) — hz sin ( 



= 0. (8.4) 



We first discuss fields along the hard-axis as they have 
the most interesting effect on the Berry phase. For 
hx,y = but hz ^ Q all four configurations that emerge 
from (jiQc in ( |3.3D are still energetically degenerate: The 
invariancc of ( |2.1| ) under 7r-rotations around the z-axis 
(which remains intact for ^ 0) implies the degeneracy 
of configurations of opposite charge but same chirality. In 
addition, states of opposite chirality and charge are also 
dege nerate since with (j){x),0{x) also —4>{x),0{x) solve 
(^yll) (with hx,y = 0). 

In the limit of large hard- axis anisotropy, ^ Ky , 
the possible ^-configurations are still (j)Qc given by (3.3) 
while 



cos 



^hz/Kz 



Inserting this into ( |8.2| ) we have 

a = a(l + hz/Kz), 



(8.5) 



(8.6) 



.1 



which demonstrates that the topological phase 
indeed be tuned by the external field. 

For arbitrary values of the ratio Kz/ Ky no analytical 
solution for the soliton structure can be found. However, 
we can convince ourselves that a is still field-dependent: 
As is verified by inserting (f>' ^ 0' = and (f> = ±7r/2 into 
( ^.4[ ), the spins far away from the soliton get pulled out 
of the easy-plane, cos0 = —hz/{2{Ky + Kz))- Thus in 
general a is different from a. 



How does the field dependence in (B.E) affect the band 
structure? Let us assume that a ~ ^aS^ is a positive 
integer, i.e. a = a — iKi/d for hz = . The dispersion 
then has the tight-binding form (6.19) of Fig. and 
consists of two coinciding chirality sheets. With increas- 
ing field hz the sheets of opposite chirality C = ± get 
separated, each shifted by Afc = |a — a|. At an external 
field /i° = Kza/2SNAd, this shift becomes Afc = 7r/2d 
and the dispersion shown in Fig. Hja is reached. Thus, a 
system with integer a can be continuously transformed, 
until it reaches half-integer behavior, and vice versa.Eil 
Note-that this behavior is periodic in the field with period 
2/i°.e3 Moreover, if the field hz{t) and thus a{t) depend 
on time, it is clear from the effective Hamiltonian ( 6.22| ) 
that da jdt plays the role of a force driving the Bloch wall 
in positive/negative a;-direction for positive/negative chi- 
rality. Note that this force has its origin in the "classical" 



part of the Berry phase, (j)cosd, and therefore can also 
be deduced from the classical Landau-Lifshitz equation 
( p. 71 ). [It is somewhat surprising that this force, as far as 
we know, has not been discussed in the literature.] 

A similar driving effect is achieved by applying an 
external field hy along the easy-axis. Indeed, insert- 
ing 4>Qc{x — X) of Eq. ( |3.3| ) into the Zeeman term 
hy J dx ainifiQc = —2hyNAQX we see that a weak mag- 
netic field acts like a (classical) force on the soliton center 
where Q is the charge of the soliton. It can be seen that 
the phase a remains unaffected by hy . Note that in anal- 
ogy to electromagnetism, where E — —A/c, hz plays the 
role of the vector potential A (albeit chirality dependent), 
while hy corresponds to the,|electric field E. Elsewhere 
we have discussed in detailElJcJ how such forces can give 
rise to Bloch oscillations of the Bloch wall — a magnetic 
analogue of the Josephson effect. Similarly, we expect a 
variety of effects for oscillating fields such as resonances 
due to the Wannier-Stark ladders and related localiza- 
tion effects. Here we just note that external fields along 
the easy or hard-axis can be used to drive the system 
through the Brillouin zone. 

Finally, we consider an external field h^ along the prop- 
agation axis. This field lifts the degeneracy between walls 
of o pposite chirality (with Q fixed), and we find from 
dJ, 2Ec = H[cl)Q^c=+] - H[cj)Q,c=-] = 'i7rQNAShx, 
which is simply the effective Zeeman splitting energjj of 
the two chirality states. From the exact solutionsE^I to 
(^.4|) we see that the phase becomes a = a + Cao{hx), 
where ao vanishes for h^ 0- Hence the relative phase 
between walls of opposite C remains 2a, independent of 
the field, and the effective Hamiltonian ( |6.22 ) becomes 



n= ^{p^aaz)^ + ViX) + eax + Ecaz. (8.7) 

Qualitatively, we see that the last term shifts the dis- 
persion sheets of opposite chirality in opposite vertical 
direction. In the free limit, V = 0, the eigenvalues are 
E± = {P+a^)/2M±[ka(ka+2MEc)/M'^- 



-El\^l\ 

Thus, the results of Sec. Vl remain basically unchanged 
for Ec^e with the level splitting at fc = becoming 
now 2y/e'^ + E^. For E g > e tunneling of the chiral- 
ity (as discussed in Sec. VII) and hence its alternation 
in the Brillouin zone will get suppressed. For instance, 
if e ~ lOmK this requires not to exceed 3 • 10^'^Oe 
(for the YIG values of Sec. VII). On the other hand, the 
field Hx provides a useful tool to enhance observability 
of the chirality switching, since it can be used to offset 
unwanted level detuning and to restore the degeneracy of 
the chirality states. 
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APPENDIX A: COHERENT STATES AND 
BERRY'S PHASE 

In this appendijp-we discuss the path integrals for co- 
herent spin stateala and, in particular, the associated 
Berry phases. We emphasize single- valuedness of spin 
states and the role of admissible gauges since this is of 
central importance for the spin parity effects discussed in 
the main text. 

A coherent state is the state of minimal uncertainty for 
spin components transverse to the spin-quantization axis. 
It is defined as the maximum eigenstate oi Sz, \S,M — 
S), rotated into the direction of the unit vector ft ~ 
(sin 9 cos 0, sin sin 0, cos 9) , 



-iS^4>-iSyd 



''\S,M^S), 



(Al) 



where S is the spin operator. By construction, the coher- 
ent state ( |Al[ ) obeys the eigenvalue equation S • r2|Jl) ~ 
S\n) and is an eigenvector o£-S^ with eigenvalue S{S+1). 
By use of Wigner's formulaEil, (Al) can be expressed as 



M=-S 
S-M 



2S 

S + M 



1/2 



9\S+M 



Sin -J \S,M). 



(A2) 



The Euler angle x to be fixed by|-t|he requirement 
that the coherent state be single value(£3 upon cj) —> (f> + 
27rn, n = 0,±1, . . .. Thus, x is only allowed to take the 
following values 



X={2n+1) 



0,±1,. 



(A3) 



For the choices n — —1 and n = we shall use the 
terms "north-" and "south-pole" gauge, respectively. Of 
course, the results obtained in either of these gauges must 
be physically equivalent. Note that this requirement of 
single valuedness has nothing to do with the transfor- 
mation properties of under active rotations by 27r 
which, of course, will always produce a factor of (— 1)^'^ 
irrespective of the choice of x- nn 

For later use we l ist a few important propertieglaEa of 
the co here nt states ( A2) in the north-pole gauge x — —<!>■ 
From (A2) it follows that coherent states are in general 
not orthogonal 



= I cos -cos - 



sin-sin-e*(*-^') 
2 2 



2S 



(A4) 



since $7 may vary continuously on the sphere while there 
are only 25* -|- 1 mutually orthogonal spin eigenstates. For 
infinitesimally separated angles, the overlap becomes 



{n'\fl) = l + iSS(P (cos6i-l). 



(A5) 



where S(j) = ((>' ~ (p. For the south pole parametriza- 
tion X = the overlap between infinitesimally separated 
states becomes 



{n'\fi) = i + iS5(t> (cose* -I- 1). 



(A6) 



Coherent states also form an overcomplete setc£l 
^^Jdn \n){fl\ = 1, where dn = dct> d{cos9). Al- 
though the states are not orthogonal, the overlap be- 
tween different states decreases for rapidly large 5' with 
increasing angle, since 



{n'\n)\ = (^^{1 + n' ■ n)^ ^ 



(A7) 



In addition we shall make use of the fact that for large 
5, we have 



(A8) 



This relation follows from the exact expressions of the 
spin matrix elements and from the fact that fluctua- 



fi'\s\fi) ^[sn + o{Vs) {^I'l^i 



tions have size 0{VS) since the overlap (A7) decreases 
as exp{-5(n' - n)^/^}. 

We derive now a path integral representation for the 
transition amplitude between two spin configurations. To 
this end, we represent the state vector of the system as a 
product of coherent states over all lattice sites = 
(8)ilfi l^i)- Following the usual procedur£§, we slice the 
interval into N identical pieces of length e and insert 
complete sets of states at each lattice site and imaginary 
time step r„ = n(3/N, 



({n,}\e-P^\{na}) - ( n n / d^^i^rn)] 

X n ({ (^n+i)}|l - eH|{f2(r„)}), (A9) 



where dfl, = ((25 + l)/4n)dn, and \{n{To)}) ^ |{r2j), 
\{fl{TN)}) = \{^b})- In the limit of large S we use ( [A^ ) 
and write 



({f2(T„+i)}|l-eW|{n(r„)}) = 

(1 - eW[5a(T„)])({0(T„+i)}|{0(T„)}), (AlO) 

where 'H[Sfli{Tn)] is the diagonal element of the Hamil- 
tonian and follows from (2.1) by substituting Si by 
S^i{T„). A 0(5^/^) correction to this diagonal el 
ment has been dropped following standard reasoning 
For large S, large deviations between coherent states 
at adjacent imaginary time steps are exponentially sup- 
pressed due to (A7). Therefore the tr ajec tories in imag- 
inary time become smooth, and from (A5) we obtain for 
the overlap between coherent states at adjacent imagi- 
nary time steps 



({f2(r„+i)}|{rj(T„)}) ~ J];{l - i5 <50,(r„)(l - cos^^,(r„))}, 



(All) 
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where 5<j)i{Tn) = 4>i{Tn+i) — 4>i{Tn)- These overlap t erm s 
are of purely kinematical origin and contribute to ( [A9| ) 
even in the absence of a Hamiltonian. It is these terms 
which are responsible for the distinct behavior of half- 
odd integral and integral spins. Passing to the time- 
continuum limit N oo we obtain, 

where 2?f2i(r) = Hn / dni{Tn) is the measure, and we 
replaced e5(j)i{T) by d<pi{T)/dT. 

In the space-continuum limit where the spin config- 
urations vary slowly over the lattice constant a the 
exchange term in H[SQi] becomes ^ ■ ^ fli ■ fli+p = 

I ^Si(^^j)^- The transition amplitude then takes 
finally the form 



({Of,}|e-^^|{fiJ) = Jvcl>V{cosl 



(A13) 



where the path integral runs over configurations that sat- 
isfy n{x,0) = fla{x), n{x,(3) = ftbix). The Euclidean 
action is given by Se = Swz + Jq^t H, where the dy- 
namics is determined by the Wess-Zumino or Berry phase 
term 



S 



/3 



d'^r / dr 4> (l - cost 



(A14) 



and the energy of the spin configuration is given by 

[ ^ I JS^a^[{V9)^ +siii^e{V(b)^] (A15) 



KyS^ sin^ sin^ (j) + K^S^ cos^ 9 



n" 



Jn' 



£if2-Aj\ 



(A17) 



over a "vector potential" Ajy^s — e^{cos9 ^ l)/sin6', 
where the upper (lower) sign corresponds to the north. 
(south)-pole parametrization. These vector potentialsEll 
are equivalent to that of a magnetic monopole of unit 
strength evaluated on the surrounding unit sphere. The 
gauge characte r of the different parametrizations of the 
coherent state (Al) becomes now apparent. If we gauge 



transform the coherent state i— > e~ \ft), where 
A = X(j), the overlap (A17) transforms according to 
Ai-^A + V^A. By the choice of the gauge one de- 
cides whether a part of the Berry phase "disappears" 
in the definition of the coherent state or whether it ap - 
pears explicitly in the path integral via the overlap ( [A17| ). 
However, in order to preserve the single- valuedness of the 
coherent states — our fundamental postulatp=n=anly gauge 
transformations exp{— iA} arc admissibleE30 that are 
single valued upon 0^0 + 2tt. Evidently, this is the 
case for A = 2S(j) (for all S) which relates the north- and 
south-pole parametrization. On the other hand, for half- 
odd integer spin this condition is violated for A = iScf) 
which relates the coherent state with x = —0 to the 
one with x = 0, but the latter is not single valued and 
thus not an admissible state. The corresponding vec- 
tor potential would be Aq = — cot 9e^ and does not 
yield the full Berry phase accumulated in a closed circuit: 
f dCl -Aq measures the area between the trajectory on the 
unit sphere and the equator while f dfl ■ Ajv^s measures 
the area between the trajectory and the north or south 
pole. For trajectories crossing the datelincEJ (this is typi- 
cally the case if spherical coordinates are chosen that are 
adapted to the symmetry of the Hamiltonian), the phase 
factor expliS" / dfl ■ Aq}, that result s fr om the "wrong" 
choice X = for the coherent state (Al), does not coin- 
cide with the Berry phase term, expliS* J dJl • Apf g}, for 
half-odd integer spins and would, e.g., lead to a wrong 
semiclassical spin-quantization. 



Eqs. (A1E)-(A16) generalize the formalism of micromag- 
netics to include quantum interference effects. Our dis- 
cussion is not restricted to the anisotropy configurations 
shown here, one could also include e.g. magnetostatic 
interactions of more general form. In the particular case 
where the spin con figu rat ion only depends on one coor- 
dinate, we recover (2^)- (2.5). 



Note that the Berry phase term ( A14) has been derived 
in the the north-pole gauge x — Ifj instead, we had 
used the south pole gauge x = we would have obtained 



wz 



S_ 
73" 



d^r / dr (/) (1 + . 



(A16) 



This gauge dependence can be traced back to the gauge 
dependence of the overlap (A5), (A6) of infinitesimally 
separated coherent states. It is instructive to express 
this overlap as a line integral 



APPENDIX B: EVALUATION OF THE 
DAMPING KERNEL 

In this Appendix we present the der ivat ion of the 
damping kernel ( 5.15 ) startin g fro m Eq. (5.8). 

In order to evaluate F in ( [5. 8] ) we first complete the 
square in the exponential. As we are working only to 
order 0{X'^ /c^), it is sufficient to shift 1^9 by p = \Q^^J 
since 

^■[Q + K]ip + J-ip={^ + p)- [g + K]{ip + p) + 0{X''). 

(Bl) 

[Q^K, are Hermitian]. Thus, Eq. ( ^.8|) can be rewritten 
as 



F[X]= / V(p det 



^dx{<t>'S-<j>l[^-p\]5{i 
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-Na0-[G+IC]v 



(B2) 



where ip = 'P + P- Eq. ( p32| ) can now be considerably sim- 
plified. First we note that J dxcji'^p cx / dx(j)QO~^(j)Q = 
due to the parity invariance of Q and (anti-) symmetry of 
4>'q {(j>o). Thus the i5-function enforces (fi to be orthogonal 
to the zero-mode. The Gaussian integrations over (f are 
then well defined and the fluctuations have effective size 
0(1/ ^/Na) ■ Rescaling ip — \/Na^P and making use of the 
identity det = exptrln, we can rewrite detiSQ / 5X) as 



exp{trln{l - ^ j cp'^, p + (5/2) J q^'^p}}, 



(B3) 



where we used J (pQ = 2/S (cf. ([3.5[)), and where the 
constant exp{trln(2/i5)} has been absorbed into the inte- 
gration measure. The second term under the logarithm 
can be neglected for large Na and the last term being 
proportional to gives rise to a pure mass renormal- 
ization. Neglecting irrelevant prefactors, Eq. (B2) thus 
becomes 



-ip-[G+K]ip 



(B4) 



where the prime on the determinant denotes omission of 
the zero mode which is enforced by the 5-function, and 
AM — 0{{Na)^) is a small mass renormalization whose 
exact value is not of interest here. In the evaluation 
of the determinant we will encounter several (ultravio- 
let) divergent terms which also have the form of a mass 
renormalization of order 0{{Na)'^)- All these renormal- 
izations will change the mass M into the experimentally 
observed "dressed" Bloch wall mass Mcs- We will thus 
drop all these renormalization terms and simply replace 
M Mcff in the action ( |57| ). 

Moreover, siijce the SG model is known to be 
renormalizableEil and since we are interested only in the 
long time (infrared) behavior there is no need here to 
carry out a systematic renormalization procedure to re- 
move the short time divergences. 

We now turn to the explicit evaluation of the deter- 
minant in (B4). We make again use of the identity 
In det = tr In and expand the logarithm: 

,-itr'iii(e[i+e-iK;]) 



^ _^^^-\tr'[g-^K-\g-^Kg-'K+o{{x/cf)) /g^N 
VdiF^ ■ 

Since JC = 0{X/c) this represents an expansion in in- 
creasing powers of X/c. The factor [det'C/]"^/^ is inde- 
pendent of X and is the partition function of spin wave 
fluctuations around the static Bloch wall. The trace in 
(B5), tr(-) = X^kC^I ■ l'*^)' evaluated in the basis of 



eigenf unctions of Q, 



e|k) =ek|k), 



(B6) 



where k = J/c^. The anisotropy gap Ky will have impor- 
tant consequences for the form of the damping kernel be- 
low. The eigenfunctions factorize into a space and (imag- 
inary) time part, |k) = |w)|A:), where {t\uj) = e*"'^/-y/9 
with Matsubara frequencies lo = 2ttv/ (3, v — 0,±1, . . .. 
Since the co ntri bution of the zero-mode (/)q is explijcijtly 
excluded in (B5), we only need the spin wave stated 



(x\k) iVfc[-ifc(5 + tanh(f)]e 



ikx 



(B7) 



where Nk = [L{1 + k'^S'^)]-^^'^ . The k values in (|B7|) are 
implicitly defined by 



kL + A{k) = 27rn, 



(B8) 



where A(fc) = 2arctan-Jj is the scattering phase shift of 
the eigenf unction (B7). 

To render the results finite in the thermodypamic limit, 
we have to subtract the vacuum fluctuationaS and thus 
renormalize, 



1 



^ detiGo+JC) 

det' (g + K.) ^ det'{g + IC)' 



(B9) 



where Qq = —kO^ — Jd^ + Ky is the operator describing 
spin waves around the anisotropy minimum in the ab- 
sence of a Bloch wall. Go has the same eigenvalues (^ ) 
as g but the space eigenfunctions are simply plane waves 
where the k values are given by fcfice — 2Tm/L rather 
than (p3^). For the results given below which only in- 
volve one summation over k, the renormalization (B9) 
then simply amounts to the replacement 



dk[p{k) - pfi-c 



— 00 

°° dkdA 



— 00 



27r dk 



(BIO) 



where p — dn/dk ~ 1^1^ density of states 

corresponding to ( |B8| ), phcc = Stt/L, and where we go 
over now to the thermodynamic limit. From the defini- 
tion of A it follows that dA/dk = -2S/{k^S^ + 1). 

With these preliminaries, we can now rewrite the low- 
est order term in (p5|) as follows. 



- Iti-g-^JC = -1 ^ - [2{k\ - id,\k) {u;\iXdr\u;) 

k 

~ {k\dl\k){u;\X'\u;)], (Bll) 



where k — J/c^. Using the eigenfunctions (B7) we obtain 

{k\-id,\k') = fc4fe, + o(L-i) {k\dl\k') = -P4fc' + 

0(L"i). Inserting the identity 1 = /dT|r)(T| we obtain 
in leading order in L, 

- ^trg^^/C = ^ E 7" / ^^{2t^^ - kX^y (B12) 
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The first term on the rhs vanishes since ek is symmetric 
in both k and lo and the second term leads to a mass 
renormahzation which d iverg es logarithmically (after the 
partial renormahzation ( BIO )). As mentioned above, this 
term is part of the dressing of the "bare" Doring mass to 
the experimentally observed value A/eff, and thus there 
is no need to remove this divergence explicitly. 

The damping due to spin waves will be exclusively due 
to the remaining terms in (B5) which shall be discussed 
next. Using the above notation, we have up to order X'^ 



1 



1 



4 ^ ' ^'^kek' 
In leading order in L we have 



(B13) 



{k\Xdrd,\k') 



kuj' 



dT e'("'-'^)^X(r). (B14) 



Thus Eq. (B13) can be rewritten in the form 



1 



4 

with 



tr {g-'KY 



7(t) 



1 



P 

dT 



E 



daX{T) X{a) -f{T ~ a), (B15) 



,'2 



(B16) 



where w| = c^(fc^ + 6 ^). With partial integrations and 
with 7(r + (3) — j{t) Eq. ( B15 ) reduces to 



itr {g-'icf = -y^r^daKir - a) [X{r) - Xia)]', 

(B17) 



where K{t) = —29^7. In ( B17 ), we have neglected a 
term 2[X{(]) — X{0)] J drXj which turns out to be small 
for typical tunneling processes. For the evaluation of 7 
and K we make use of the exact relation 



00 

E 



cosh(ti;(|i 



D) 



sinh 



(BI8) 



where w„ = 2?™//? and the rhs is understood to be pe- 
riodically extended beyond |t| < (3/2. With (B16) and 

-2dh 



(B18) we finally obtain for K 



2Dl(r) . (B19) 



Note that as a consequence of the relevant coupling be- 
tween the system X and bath which is quadratic in 
the bath coordinates (/?, K is proportional to Z?J rather 
than as in the usual Caldeira-Leggett theory. For 
low temperatures, the damping kernel (B19) reduces to 
K{T) = -2Y.uk^ule-^-^\^\. 
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FIG. 2. a) Bloch wall configuration with Q = 1, C = — 1 
in a thin long slab centered at the pinning site X — 0; h) 
Periodic pinning potential V for the wall center X. 



FIG. 3. Dispersion of a soliton in a weak periodic po- 
tential, a) For a — NASd/a integer the dispersion resembles 
that of a Bloch electron and the gaps at ±n/d are due to the 
periodic pinning potential (3.3). b) For a half-odd integer 



the Brillouin zone is halved and two subsequent band minima 
have opposite chirality. Band gaps 2e arise due to tunneling 
between the chiralities as de scribe d by ( 6.22 ), E+ and E- are 
the dispersions as given by (3.23). 



FIG. 4. Soliton dispersion in the tight binding limit, a) 
For a = NASd/a integer a standard tight binding dispersion 
results, b) For a half-odd integer, the Brillouin zone (and 
bandwidth) is halved and two subsequent band minima be- 
long to opposite chiralities. A gap 2e develops if the two 
chiralities of the soliton are connected by tunneling. 



FIG. 1. The Berry phase factor for one single spin S, 
exp{iS f dcj>{l — cos 9)} = e'^^, where A is the area on the 
unit sphere enclosed by the trajectory C traced out by S. 
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